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ABSTRACT 

We report detection of thermal emission from the exoplanet WASP-19b at 3.6, 4.5, 5.8 
and 8.0 u\m. We used the InfraRed Array Camera on the Spitzer Space Telescope to 
observe two occultations of WASP-19b by its host star. We combine our new detections 
with previous measurements of WASP-19b's emission at 1.6 and 2.09 u.m to construct 
a spectral energy distribution of the planet's dayside atmosphere. By comparing this 
with model-atmosphere spectra, we find that the dayside atmosphere of WASP- 19b 
lacks a strong temperature inversion. As WA SP-19 is an active star (log B!„y = — 4 .50± 
0.03), this finding supports the hypothesis of lKnutson. Howard fc Isaacson! (|2010f ) that 
inversions are suppressed in hot Jupiters orbiting active stars. The available data are 
unable to differentiate between a carbon-rich and an oxygen-rich atmosphere. 

Key words: methods: data analysis - techniques: photometric - occultations - plan- 
ets and satellites: atmospheres - planets and satellites: individual: WASP-19b - stars: 
individual: WASP-19. 



1 INTRODUCTION 

By observing the occultation of an exoplanet by its host star, 
we can measure the emergent flux from the planet's day- 
side atmosphere. Such measurements are challenging due to 
the low planet-to-star flux ratio (for the best cases, typ- 
ically a few tenths of one per cent in the near infrared) 
and s ources of noise both instrumental and stellar in origin 
(e.g.: iKnutson et"aI1l2008l : [s~mith et al.ll20riT). To date, most 
such measurements (inclu ding the first, Charbonneau et al.l 
120051 : iDeming et al.|[2005h have been made with the Spitzer 
Space Telescope, though ground-base d facilities are now 
an essential complement ( e.g.: Ide Mooii &: Snellen! 120091 : 
ISing fc L6pez-Moralesll2009l ). 
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With a single occultation we can measure the cor- 
responding brightness temperature an d determine the ec- 
centr icity of the planet's orbit (e.g. ICharbonneau et alj 
■ 2005), which is necessary for the accurate determination 
of the stellar and pla netary radii in transiting systems 
ijAnderson et al.ll2011al ) and is important for studies of the 
formation and tidal inflat i on of short-period, giant plan- 
ets ( e.g.: iNagasawa fc Idal l201ll ; llbgui. Burrows fc Sp iegel 
120101 ) 



With photometric measurements at various wave- 
lengths, we can construct a spectral energy distribution of 
the planet's dayside atmosphere. From this we can infer 
properties such as the planetary albedo, the dayside en- 
ergy budget and th e efficiency of dayside-to-nightside en ergy 
redistribution (e.g. iBarman. Hauschildt fc Allardll2005l ). As 
the atmospheric depth probed depends on the molecular 
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opacity in the observation band, a spectrum can describe 
an atmosphere's vertical temperature structure. The obser- 
vation that some planet atmospheres exhibit strong tem- 
pera ture inversions, or stratospheres, a nd others do not 
(see iKnutson. Howard fc Isaacson! 120101 . for a summary) 
led to a suggestion that inversions are present in atmo- 
spheres hot enough to maintain high-opacit y absorbers in 
the g as phase in the upper atmosphere (e.g. iFortnev et alj 
2008). However, this was challenged by recent contrary re- 



sults jMadralekeUalJ|2^^ and theory 

(e.g. Spiegel, Silvcrio & Burrows 2009). As an alternative, 
IKnutson. Howard fc Isaacson! (|2010h suggested that those 
planets orbiting chromospherically active stars lack inver- 
sions because the associated high UV flux destroys the high- 
opacity, high-altitude compounds that would otherwise in- 
duce inversions. This hypothesis is based on a small sample 
and further measurements across a wider parameter space 
are vital to test it. Specifically, there is a paucity of mea- 
surements for planets orbiting active stars. We suggest com- 
position could be a key factor, with low-metallicity planets 
lacking high concentrations of the high-opacity absorbers, 
such as TiO and sulphur, thought to be responsible for in- 
versions (e.g. lFortnev et alj 12008 : IZahnle et alj|2009l ). 

As the spectral coverage increases, so too does the 
information about a n atm osphere that we can discern. 
iMadhusudhan et alj (|201 lal) used seven measurements of 
the e mission of WASP- 12b (|Croll et al.ll201ll : ICampo et all 
2011), to show that the dayside atmosphere is the first 
known to be carbon-dominated. Thus we are entering an 
era in which we can make statistically-sound inferences 
about the composition of e xoplanet atmospheres. Further, 
IMadhusudhan et alj (|2011af) demonstrated that the planet 
lacks a prominent thermal inversion and has very inefficient 
day-night energy circulation. 

In this paper, we present Spitzer measurements of 
WASP-19b's dayside thermal emission at 3.6, 4.5, 5.8 and 
8.0 um. Discovered by the Wide Ang le Search for Planets 
l|PoUacco et al.l 2006: He bb et al.ll201oh . the shortest-period 
hot Jupiter, WASP-19b, is a 1.17-Mj u p planet in a near - 
circular, 19-hr orbit around a G8V star. lHellier et all (|201ll ) 
showed the planet's orbital axis to be aligned with the spin 
axis of its host star. The planet's ther mal emission was pre- 
viously measured a t 1.6 and 2.09 [im (jAnderson et al.ll2010l : 
iGibson et afll201Ch . thus we bring the measurement tally to 
six bands. We use all six thermal emission measurements 
to characterise the planet's atmosphere and refine the sys- 
tem parameters by combining these data with pre-existing 
photometry of the transit and radial- velocity data. 



2 NEW OBSERVATIONS 

We observed two occultations of the planet WASP- 19b by 
its host star WASP- 19 (2MASS 1559509 5-2803422, K s = 
10.22) with Spitzer jWerner et all 120041 ') during UT 2009 
January 29 and UT 2009 March 22. On each date, we em - 
ployed the Infrared Array Camera (IRAC, iFazio et al 1120041 ) 
in full array mode (256 x 256 pixels, 1.2 " pix -1 ). During 
the first occultation, we measured the WASP-19 system si- 
multaneously in the 4.5- and 8.0-u.m channels (respectively, 
channels 2 and 4) for a duration of 3.2 h. We measured the 
second occultation simultaneously in the 3.6 and 5.8 urn 



channels (respectively, channels 1 and 3) for a duration of 
3.0 h. 

Prior to the first occultation, we used the emission neb- 
ula NGC 7538 in Cepheus to 'pre-flash' the target position 
on the detector arrays for 0.5 hr. This was intended to reduce 
or remove the known illumination-history de pendence of the 
gain response of the 8.0-u.m detector (e.g. IKnutson et alj 
120081 . and references therein). The 5.8-um detector is known 
to suffer a similar issue, but pre-flashing this array was not 
permitted due to the detrimental effect it is known to have 
on the array. Instead, we attempted to stabilize the array 
by observing the target for an extra hour prior to the oc- 
cultation. As this 'pre-stare' was performed as a separate 
observation request to the occultation observation, the tar- 
get was reacquired between the two. As such, the pre-stare 
observation must be treated as a separate dataset with its 
own systematics and so are of little use in determining the 
occultation depth. 

Using an effective integration time of 10.4 s, we ob- 
tained 876 and 840 images, respectively, for the first and 
second occultations. In each dataset, we see a small, peri- 
odic, pointing wobble (P ~ 1 hr), thought to be caused by 
the thermal cycling of an on-board battery heateiQ. There 
is also a very small drift of the target position over the span 
of each dataset. 

We used the images calibrated by the standard Spitzer 
pipeline (version S18.7.0) and delivered to the community as 
Basic Calibrated Data (B CD). Our method is ess entially the 
same as we presented in I Anderson et all (|201ld ). to which 
we refer the reader for further information. For each im- 
age we converted flux from MJy sr _1 to electrons and then 
used iraf to perform aperture photometry for WASP-19, 
using circular apertures with a range of radii: 1.5-6 pixels 
for the 3.6 um and 4.5 [im data and 1-5 pixels for the 5.8 
Um and 8.0 um data. The apertures were centred by fit- 
ting a Gaussian profile on the target. The sky background 
was measured in an annulus extending from 8 to 12 pixels 
from the aperture centre, and was subtracted from the flux 
measured within the on-source apertures. We estimated the 
photometric uncertainty as the quadrature addition of the 
uncertainty in the sky background (estimated as the stan- 
dard deviation of the flux in the sky annulus) in the on- 
source aperture, the read-out noise, and the Poisson noise 
of the total background-subtracted counts within the on- 
source aperture. We calculated the mid-exposure times in 
the HJD (UTC) time system from the MHJD.OBS header 
values, which are the start times of the DCEs (Data Col- 
lective Events), by adding half of a DCE duration (FRAM- 
TIME). 

The choice of aperture radius for each dataset was a 
compromise between maximising the signal-to-noise of the 
measurements and, from fits to all available data (see Section 
3), minimising the residual scatter of the lightcurve. Each 
consideration suggested very similar optimal apeture radii 
and we adopted 2.7 pix for the 3.6- and 4.5-(im data and 
2.5 pix for the 5.8- and 8.0-nm data. For each dataset we 
found that the variation in the fitted occultation depth (see 
Section 3) was much smaller than 1 a over a wide range of 
aperture radii. 



http: / /ssc. spitzer. caltech.edu/warmmission/news/21oct2010memo.pdf 
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Table 1. Number of points rejected per dataset per criterion 



Datasct 


Flux 


x-pos 


y-pos 


Total (%) 


2009 Mar 22 / 3.6 (im 


8 


1 


3 


11 (1.3) 


2009 Jan 29 / 4.5 urn 


11 


11 


11 


51 (5.8) 


2009 Mar 22 / 5.8 [im 


11 


9 


7 


21 (2.5) 


2009 Jan 29 / 8.0 urn 


4 


8 


9 


17 (1.9) 



Some groups choose to reject a portion of data at the 
beginning of each observation, citing as justification either 
the settling of the spacecraft (not seen in our data) or an 
improvement in the fit; we found no reason to do this. We 
rejected any flux measurement that was discrepant with the 
median of its 20 neighbors (a window width of 4.4 min) by 
more than four times its theoretical error bar. We also per- 
formed a rejection on target position. For each image and 
for the x and y detector coordinates separately, we com- 
puted the difference between the fitted target position and 
the median of its 20 neighbors. For each dataset, we then 
calculated the standard deviation, a, of these median differ- 
ences and rejected any points discrepant by more than 4cr. 
The numbers of points rejected on flux and target position 
for each dataset are displayed in Table [T] According to the 
IRAC handbook, each IRAC array receives approximately 
1.5 solar-proton and cosmic-ray hits per second, with ~2 
pixels per hit affected in channels 1 and 2, and ~6 pixels 
per hit affected in channels 3 and 4, and the cosmic ray flux 
varies randomly by up to a factor of a few over time scales 
of minutes. Thus, the average probability per exposure that 
pixels within the stellar aperture will be affected by a cos- 
mic ray hit is 1.3 per cent for channels 1 and 2 and 3.2 per 
cent for channels 3 and 4, which is in good agreement with 
the portion of frames that we rejected. These probabilities 
are likely to be underestimates as we calculated them using 
partial pixels and neglecting the effect of hits within the sky 
annuli. For an unknown reason, a greater portion of channel 
2 images were rejected due to jumps of the target position, 
particularly in the x direction. The post-rejection data are 
displayed raw and binned in the first and second panels re- 
spectively of Figure [1] 



3 DATA ANALYSIS 
3.1 Data and model 

We performed a global determination of the system pa- 
rameters incorporating: our new Spitzer occultation pho- 
tometry; the HAWK-I Jz~-band an d 2.09-uMn occu l tation 
ligh tcurves obtain e d, res pectively, bv lAnderson et all (|2010l ) 
and iGibson et all l|201Ch ; the 34 CORALIE radia l-velocity 
(RV) measurements listed in iHebb et al.l (|2010l ); the 36 
HARPS RVs, obtained through a transi t, and the 3 
CORALIE RVs given in IHellier et al] |201ll); the L COGT 
FTS z-band transit lightcurve from IHebb et all (2010); 
and the ESQ NTT r-band transit lightcurve presented in 
IHellier et aI1l|201ll) . We did not include the three seasons of 
WASP survey photometry presented in IHebb et al.l (2010). 
Rather, we placed a Bayesian Gaussian prior on the epoch 
of mid-transit, T c , using the epoch given in IHellier et alj 



i|201ll ): T c = 2455168.96801 ± 0.00009 HJD. Thus our anal- 
yses completed quicker and the shape of the transit was 
determined using only the high-S/N photometry. This can 
be preferable as photometry from surveys such as WASP is 
prone to dilution and, depending on which detrending al- 
gorithm is used, the transit depth can be suppressed. We 
decorrelated each transit lightcurve with a linear function 
of time. The H AWK-I JJ-band data w ere partitioned and 
detrended as in lAnderson et~ai] l|2010l ) and the HAWK-I 
2.09-um ligh tcurve was decorrlela ted with a linear function 
of time as in IGibson et ail l|2010l ). These data were used as 
input into an adaptive Markov-chain Monte Carlo (MCMC 
algorithm (ICollier Cameron et al.l2007l ; Pollacco et al .1200; 



lEnoch et al.l201oh ; see Anderson et al.l ( 201 lbl) for a descrip- 
tion of the current version of our code. Such an analysis, in- 
corporating all available data, is necessary to take account 
of the cross-dependancy of system parameters and to make 
a reliable assessment of their uncertainties. We partitioned 
the RV data by spectrograph so as to allow for an instru- 
mental offset and for a potential specific stellar activity level 
during the short-baseline HARPS observations. 

The MCMC proposal parameters we used are: T c , 
P , (Rpi /R*) 2 , Ti4, 6, Ki, T off , [Fe/H] , </k~ cosoj, -yesincj, 
\/v sin i cos A, Vv sin i sin A, AFi.e, AF2.09, AF3.6, AF4.5, 
AF5.8, A_Fg.o and i ff • See Section 13.21 for a definition of 
t g and Table [3] for definitions of the other parameters. At 
each step in the MCMC procedure, each proposal parame- 
ter is perturbed from its previous value by a small, random 
amount. From the proposal parameters, model light and RV 
curves are generated and \ 2 i s calculated from their com- 
parison with the data. A step is accepted if \ 2 ( our merit 
function) is lower than for the previous step, and a step 
with higher % 2 ls accepted with probability exp(— A^ 2 /2). 
In this way, the parameter space around the optimum so- 
lution is thoroughly explored. The value and uncertainty 
for each parameter are respectively taken as the median 
and central 68.3 per cent confidence interval of the param- 
eter's ijnarginalised posterior probability distribution (e.g. 
lFordll2006l ). 

3.2 Spitzer data 

IRAC uses an InSb detector to observe at 3.6 and 4.5 u.m, 
and the measured flux exhibits a strong correlation with the 
position of the target star on the array. This effect is due to 
the inhomogeneous intr a-pixel sensitivity of t he detector and 



. well-documented (e.g. Knutson et al 



therein). Following I Charbonneau et al 



2008, and references 



(2008) we modelled 



this effect as a quadratic function of the sub-pixel position of 
the PSF centre, with the ad dition of a cross-term to permit 
rotation l|Desert et al.ll2009P ) and a linear term in time: 



df = ao+a x dx+a y dy+a xy dxdy+a xx dx 2 +a yy dy 2 +a t dt (1) 

where df = f — f is the stellar flux relative to its weighted 
mean, dx = x — x and dy — y — y are the coordinates of the 
PSF centre relative to their weighted means, dt is the time 
elapsed since the first observation, and ao, a x , a y , a xy , a xx , 
a yy and at are coefficients. 

IRAC uses a SiAs detector to observe at 5.8 and 8.0 
urn, and its response is thought to be homogeneous, though 
another systematic affects the photometry. This effect is 
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Figure 1. In each of the above three plots, from top to bottom, the data were taken at 3.6, 4.5, 5.8 and 8.0 um. Relative flux offsets 
were applied to datasets for clarity. Left: Raw Spitzer data with the best-fitting trend and occultation models superimposed. Middle: 
The same data binned in phase (A(f> = 0.003 ~ 3.4 min) with the best-fitting trend models superimposed. Right: The binned data 
after dividing by the best-fitting trend models, and with the best-fitting occultation models superimposed. The error bar on each binned 
measurement in the panels in the middle and on the right is the standard deviation of the points within the bin. 



Table 2. Trend model parameters and coefficients 





3.6 um 


4.5 urn 


5.8 um 


8.0 urn 


/ 


151194.55 


80580.23 


12693.57 


18991.68 


X 


31.93 


24.50 


25.99 


25.44 


9 


24.91 


25.98 


25.16 


23.63 


ao 


-24 ±16 


89.7 ±5.0 


-95.2t^ 


1569±1« 


O-x 


6566 ± 75 


-433 ± 36 






a y 


9278 ± 49 


4062 ± 32 






( ' .r y 


6040 ± 3280 


-7860 ± 4150 






&XX 


-6330 ± 1230 


23000 ± 880 






"yy 


-22500 ± 1980 


30360 ± 3810 






at 


1050 ± 270 


-2203 ± 57 






ai 






-52± u 




a.1 






-s-9±*:i 


78±J| 


toff /min 






10.6t£ 4 


15.3+^V 



ni- 
an 



known as the 'ramp' because it causes the gain to 
crease asymptotically over time for every pixel, with 
amplitude dependin g on a pixel's illumination history (e.g. 
iKnutson et al. 20081. and refe rences therein). Again follow- 
ing ICharbonneau et alj (|200ct ). we modelled this ramp as a 
quadratic function of ln(dt): 



where t g is a proposal parameter restricted to positive val- 
ues. To prevent t n from drifting to values greater than an 
hour or so, we place on it a Gaussian prior by adding a 
Bayesian penalty to our merit function (x 2 ) : 



Bp t<,tf = *off /cr t 2 off 



(3) 



where at 



15 min. 



df = a + ai ln(dt + t off ) + a 2 (ln(dt + t oS )) 2 



(2) 



A steep ramp is evident in the 8.0-u.m lightcurve (Fig- 
ure [1] middle panel). Due to the large distance on the sky 
between the target and the pre-flash source, there was an 
11-minute gap between the end of the pre-flash observations 
and the start of the target observations. It may be that the 
detector traps de-populated during this time, giving rise to 
the observed ramp that is more typical of observations with- 
out pre-flash. 

In addition to Equation we tried trend functions with 
a variety of time dependencies: no time dependency; a linear- 
logarithmic time dependency (equivalent to setting a 2 = 
in Equation [2J) ; a linear time dependency; and a quadratic 
time dependency. Each of these functions result in depths 
consistent within 1-a with the depths obtained using Equa- 
tion [5] For this reason and because Equation [5] has been 
shown to accurately describe the ramp in higher cadence, 
longer baseline datasets obtained with the SiAs detectors 
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fe.g. lKnutson et al.l l2009l). we adopt Equation[2]as our trend 

model for channels 3 and 4. 

W e used singular value decomposition (|Press et al.l 
Il992ft to determine the trend model coefficients by linear 
least-squares minimization at each MCMC step, subsequent 
to division of the data by the eclipse model. The best-fitting 
trend models are superimposed on the binned photometry 
in the middle panel of Figure [1] Table [2] gives the best- 
fitting values for the trend model parameters and coefficients 
(Equations [1] and [2} , together with their 1-er uncertainties. 

3.3 Photometric and RV noise 

We scaled the formal photometric error bars so as to obtain 
a reduced % 2 °f un ity, applying one scale factor per dataset. 
The aim was to properly weight each dataset in the simulta- 
neous MCMC analysis and to obtain realistic uncertainties. 
The error bars of the FTS and the NTT photometry were 
multiplied, respectively, by 0.71 and 1.15. The scale factors 
for the error bars of the occultation photometry from 1RAC 
channels 1, 2, 3 and 4 were, respectively, 1.02, 0.99, 1.16 and 
1.04. Importantly, the error bars of the occultation photom- 
etry were not scaled w hen deciding wh i ch tre nd models or 
aperture radii to use. In lAnderson et all lj20ld ). the HAWK- 
I H-b&nd occultation data were split into eleven lightcurves 
according to offset and telescope repointing. Within a global 
analysis, each lightcurve was detrended individually and the 
error bars of each lightcurve was rescaled by its own fac- 
tor. As the number of datapoints in each lightcurve is small 
(8-12 in the 6 lightcurves obtained prior to repointing and 
40-44 in the 5 lightcurves obtained post repointing), we here 
opted to rescale the error bars of all eleven lightcurves by the 
same factor, which was was 1.14, though we did still detrend 
each lightcurve separat ely. The erro r bars of the 2.09-u.m 
light curve presented in |Gibson et all l|2010h were scaled by 
6.15. iGibson et all (|201oF T also found their uncertainties re- 
quired a large scaling factor (6.17). They attributed this to 
variations in the inter-pixel sensitivity that likely resulted 
from their random dithering (radius = 30") of the point- 
ing between observations. Thoug h our //-band observa tions 
were also made using HAWK-I l|Anderson et al.ll201Ch . the 
corresponding scale factor was much closer to unity. This 
is because we dithered over a fixed pattern of six offsets 
and did not employ random jitter. We could thus produce 
one lightcurve per offset position and model each of their 
systematics independently. Current best practice is to avoid 
offsetting at all. 

We assessed the presence of correlated noise in the oc- 
cultation and transit data by plotting the rms of their binned 
residuals (Figure[2]). The 4.5-u.m occultation lightcurve and 
the FTS transit lightcurve each display a small amount of 
correlated noise on timescales of 5 min and longer. 

To obtain a reduced spectroscopic-^ 2 of unity and to 
balance the different datasets in the MCMC, we added in 
quadrature a jitter of 14.1 ms _1 to the uncertainties of the 
CORALIE RVs and 6.9 ms " 1 to the un c ertain ties of the 
HARPS RVs, as was done in lHellier etafl l|201lf ). 



3.4 Time systems and light travel time 

The Spitzer, HAWK-I, FTS and NTT photometry are in the 
HJD (UTC) time system. The CORALIE and HARPS RVs 



o 
CO 




0.01 
0.001 
0.0001 

0.01 
0.001 
0.0001 



Bin width / minutes 

Figure 2. RMS of the binned residuals for, from top to bottom, 
the pre-existing HAWK-I occultation photometry at 1.6 urn and 
2.1 um, the new Spitzer occultation photometry at 3.6, 4.5, 5.8 
and 8.0 um, and the pre-existing FTS and NTT transit photome- 
try. The solid black lines, which are the rms of the unbinned data 
scaled by the square root of the number of points in each bin, 
show the white-noise expectation. The ranges of bin widths are 
appropriate for the datasets' cadences and durations. 



are in the BJD (UTC) time system. The difference between 
BJD and HJD is less than 4 s and so is negligible for our 
purposes, and the timing information mostly comes from the 
photometry. Leap second adjustments are made to the UTC 
system to keep it close to mean solar time, so one should 
really use Terrestrial Time. However, our observations span 
a short baseline (2008-2010), during which there was only 
one leap second adjustment, so our choice to use the UTC 
system has no impact. 

The occultation of WASP- 19b occurs farther away from 
us than its transit does, so we made a first order correction 
for the light travel time. We calculated the light travel time 
between the beginning of occultation ingress and the begin- 
ning of transit ingress to be 15 s. We subtracted this from 
the mid-exposure times of the Spitzer and HAWK-I occul- 
tation photometry. For comparison, we measure the time of 
mid-occultation to a precision of 35 s. 
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3.5 Results 

Table [3] shows the median values and the 1-a uncertainties 
of the fitted proposal parameters and derived parameters 
from our final MCMC analysis. Figure [T] shows the best- 
fitting trend and occultation models together with the raw 
and detrended Spitzer data. Table [2] gives the best-fitting 
values for the parameters of the trend models, together with 
their 1-a uncertainties. Figure [3] displays all the photometry 
and RVs used in the MCMC analysis, with the best-fitting 
ecli pse and radial- ve l ocity models superimposed. 

lAnderson et all (|201dh measured an occultation depth 
of 0.259 ± 0.045 per cent from their //-band lightcurve. In 
analysing the same data, we obtained a similar depth of 
0.276 ±0.044 per cent. The minor difference is probably due 
to some combination of the slightly different manner in how 
the lightcurves' error bars were rescaled (see Section I3.3[) 
and the slight difference in the occultation ephemeris (the 
additional radial-velocity and occultation data result in a 
mid-occultation time at the ti me of the -ff-band obser vations 
~5 min earlier than found by lAnderson et all |2010l )). The 
occultation depth of 0.366 ± 0. 067 per cent that we derived 
from the 2.09-|J.m lightcurve of iGibson et~"aH (|2010l ) is near- 
identical to the depth that they obtained (0.366 ± 0.072 per 
cent). 

We calculated the brightness temperatures that corre- 
spond to the measured occultation depths and present these 
in Table [3] To calculate these, we defined the product of the 
planet-to-star area ratio and the ratio of the band-integrated 
planet-to-star flux densities, corrected for the wavelength- 
dependency of the transmission^, to be equal to the mea- 
sured occultation depth (e.g., Charbonneau et al. 2005). We 
used a model spectrum of a star with the T e g, logfl an d 
[Fe/H] values of TabIel3l l|Hauschildt. Allard fc Baronlll999h . 
normalised to reproduce the integrated flux of a black body 
with T e ff = 5475 K. The uncertainties in the brightness tem- 
peratures only take into account the uncertainties in the 
measured occultation depths. 

To obtain reliable determinations of the occultation 
depths and orbital eccentricity, it is important that the 
time of mid-transit is known with accuracy and preci- 
sion at the epochs the occultation data are obtained. The 
WASP photometry, which we used to place a prior on 
the transit ephemeris, span a two-year baseline of 2006 
May to 2008 May and the FTS transit lightcurve was ob- 
tained in 2008 December. All occultation data were obtained 
soon after: during 2009 January to April. The NTT tran- 
sit lightcurve, obtained in 2010 February, ensured a reli- 
able transit ephemeris at the occultation epochs. We note, 
though, that the difference between the transit ephemeris 
presented herein a nd that presented in the discovery paper 
|Hebb et al.ll2010l . i.e. without the NTT lightcurve), propa- 
gated to the occultation epochs, is a mere ~20 s. The ac- 



curacy of the discovery-paper ephemeris is due to the long 
baseline of the WASP photometry and the high quality of 
the FTS lightcurve, which was obtained only months before 
the occultation data. 



3.6 Stellar activity 

iHebb et al.l l|2010l ) reported a rotational modulation of the 
WASP lightcurves with a period of 10.5 ± 0.2 days and an 
amplitude of a few mmag. This indiciated that WASP- 19 is 
an active star, with the sinusoidal modulation being induced 
by a non-axisymmetric distribution of starspots. 

We determine the log7?HK activity index of WASP- 19 
by measuring the weak emission in the cores of the Ca 
ii H+K lines (iNoyes, Weiss fc Vaughanl [l984l ; ISantos et al] 
l200d; iBoisse et ; aUl2009h . The 36 HARPS spectra presented 



Hellier et al.l (|201lf ) had SNR in the range 14-38. We se- 



2 For the HAWK-I measurements, the transmission of 
the atmosphere, telescope, instrument, and detector were 
accounted for by using the transmisson curve obtained 
from http://www.eso.org/observing/etc/ For the Spitzer 
measurements, the telescope throughput and detector 
quantum efficiency were accounted for by using the the 
full array average spectral response curves available at 
http : //irsa. ipac . caltech . edu/data/SPITZER/docs/irac/calibrat 



lected the 12 spectra with SNR>19 per pixel at 550 nm, as 
the activity level tends to be systematically under- or over- 
estimated for spectra with low SNR. By assuming B — V 
= 0.570, we infer log R^k = -4.50 ± 0.03, which is the 
weighted mean and standard deviation of the values de- 
termined from individual spectra; we used the SNR as the 
weighting factor. T his is similar to the value of log R'kk = 
—4.66 measured bv lKnutson. Howard fc Isaacson! (|2010h . It 
is difficult to judge the level at which t he two values agree 
as iKnutson. Howard fc Isaacson! ([2010) do not provide an 
uncertainty estimate and our uncertainty value is likely to 
be an underestimate. 

As we know the true stellar rotation period to be 
10.5 ± 0.2 days from rotational modulation, we can use our 
log-RnK value to test the act ivity-rotation calibration of 
iMamaiek fc Hillenbrand! (2008). The calibration suggests a 
stellar-rotation period of P Iot = 12.3 ± 1.5 d, which is con- 
sistent within errors. 

We considered whether stellar variability could have af- 
fected our measured occultation depths. One potential issue 
is that that the stellar brightness may have varied signif- 
icantly during one or more of the observations. However, 
with observation durations of ~3 hr and a stellar rotation 
period of 10.5 d, the visible portion of the stellar surface will 
have changed little during any one observation. To first or- 
der, the resulting small impact on the occultation lightcurves 
can be modelled as a linear trend, which will be handled 
by the trend functions. Another concern is that the stel- 
lar brightness may have changed significantly between the 
non-simultaneous occultation observations. For example, the 
3.6-(i.m data were obtained two months after the 4.5-(xm 
data and it is the relative measurements at these two wave- 
lengths that are the prime diagnostic for the presence of 
an atmospheric temperature inversion. Assuming a constant 
planet brightness, the stellar brightness would need to have 
changed by ~5 per cent to have changed the occultation 
depth by 1 a and the amplitude of the modulation of the 
WASP lightcurves (a few mmag) shows that this is very un- 
likely. Thus our derived eclipse depths, and the conclusions 
on which they depend, are insensitive to the variability of 
ion44ASf)'-dpjectralresponse/ 
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Table 3. System parameters of WASP-19 

Parameter Symbol Value Unit 



Orbital period 

Epoch of mid-transit (HJD, UTC) 

Transit duration (from first to fourth contact) 

Duration of transit ingress i=a duration of transit 

Planet-to-star area ratio 

Impact parameter 

Orbital inclination 

Semi-amplitude of the stellar reflex velocity 

Centre-of-mass velocity 

Offset between HARPS and CORALIE 

Orbital eccentricity 
Argument of periastron 



Phase of mid-occultation, having accounted for light travel time 
Occultation duration 

Duration of occultation ingress sa duration of occultation egress 



P 

T c 

T14 

T12 ~ T34 

(Rpi/R*) 2 
b 



Ki 

7rvl 

THARPS — CORALIE 



ecosoj 
e sin £j 

^mid — occ. 



0.78883951 ±0.00000032 
2455183.16711 ± 0.000068 
0.06549 ± 0.00035 
0.01346 ± 0.00052 
0.02050 ± 0.00024 
0.656 ±0.015 
79.42 ± 0.39 

256.9 ±2.7 
20 787.26 ±0.23 
19.32 ± 0.82 



0.0019 1 



1-0.0085 
-0.0015 

< 0.027 (3 a) 



T 58 

T56 ~ T 78 



75 -162 

00010+ 00081 

U.UUU1U_ 00069 

0007+ 0062 
u - uuu ' -0.0016 

50006+ - 00052 

U - 0UUUD -0. 00044 

06564 +<u " ,ul:1 
u.uuom_ .00040 

0.01358 ± 0.00054 



Relative planet- to-star flux at 1.6 |J.m 


AF1.6 


0.276 ± 0.044 


per 


cent 


Relative planet- to-star flux at 2.09 (im 


AF2.09 


0.366 ± 0.067 


per 


cent 


Relative planet- to-star flux at 3.6 (im 


AF3.6 


0.483 ±0.025 


per 


cent 


Relative planet-to-star flux at 4.5 [im 


AF4.5 


0.572 ± 0.030 


per 


cent 


Relative planet- to-star flux at 5.8 [im 


AF5.8 


0.65 ±0.11 


per 


cent 


Relative planet- to-star flux at 8.0 (im 


AFg.o 


0.73 ±0.12 


per 


cent 


Planet brightness temperature^ at 1.6 \im 


1\>,1.6 


2750 ± 130 


K 




Planet brightness temperature^ at 2.09 |xm 


^b,2.09 


2670 ±170 


K 




Planet brightness temperature^ at 3.6 [im 


7b,3.6 


2346 ± 57 


K 




Planet brightness temperature^ at 4.5 um 


Tb,4.5 


2273 ± 64 


K 




Planet brightness temperaturet at 5.8 [im 


Tb,5.8 


2260 ± 230 


K 




Planet brightness temperaturet at 8.0 [im 


Tb,8.0 


2260 ± 250 


K 





Sky-projected stellar rotation velocity 


v sini 


4.63 ±0.27 


km s 


Sky-projected angle between stellar spin and planetary orbit axes 


A 


4.1 ±5.2 





Star mass 


M* 


0.969 ± 0.023 


M Q 


Star radius 


R* 


0.993 ±0.018 


R@ 


Star density 


P* 


0.990 ± 0.043 


PQ 


Star surface gravity 


logfl* 


4.430 ± 0.012 


(cgs) 


Star effective temperature 




5475 ± 98 


K 


Star metallicity 


[Fe/H] 


0.02 ±0.09 


(dcx) 


Planet mass 


M pl 


1.165 ±0.023 


Mj up 


Planet radius 


R P i 


1.383 ±0.031 


Rjup 


Planet density 


P P i 


0.440 ± 0.026 


Pjup 


Planet surface gravity 


loggp 


3.133 ±0.017 


(cgs) 


Semi-major axis 


a 


0.01653 ±0.00013 


AU 


Planet equilibrium temperaturet (full redistribution) 


Tp,A=0,f=l 


2045 ± 41 


K 


Planet equilibrium temperaturet (day side redistribution) 


Tp,A=0,f=2 


2432 ± 49 


K 


Planet equilibrium temperaturet (instant re-radiation) 


?P,A=0,f=8/3 


2614 ±52 


K 



' We modelled both star and planet as black bodies and took account of only the occultation depth uncertainty, which dominates. 
* Tp : A=o,f = / 4 Teff\/ -j 4, where / is the redistribution factor, with / = 1 for full redistribution, / = 2 for dayside redistribution 
and / = 8/3 for instant re-radiation (Cowan & Agol 2011). We assumed the planet albedo to be zero, A = 0. 
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0.93 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.4 0.45 0.5 0.55 0.6 0.6 0.8 1 1.2 1.4 

Orbital Phase Orbital Phase 

Figure 3. The results of our global analysis, which combines our new Spitzer occultation photometry with pre-existing transit photom- 
etry, ground- based occultation photometry and radial velocities. The model s generated from the best-fitting pa rameter values of T able [3] 
are overplotted. Left: From top to bottom, HAWK-I occultations at 1.6 |^m jAnderson et al.ll201(t ) and 2.09 \xm llGibson et aj201Ch . The 
error bar on ea ch binned measur ement is the standard deviation of the points within the bin. Top-right: z-band transit lightcurve taken 
with FTS ftop; lHebb et alj|201ol ) and r-band transit lightcurve taken with NTT (bottom; iHellier et alfeoilf) . Arbitrary offsets have been 
applied to the photometry plot s for display purpo ses and Bottom-right: Spectroscopic orbit and transit illustrated by CORALIE and 
HARPS data jHebb et alj |2010; Hcllicr et al. 201l|). The measured systemic velocities (Tabled of each dataset have been subtracted. 



4 DISCUSSION 



4.1 Atmosphere model 



We interpret our observations of the hot Jupiter WASP- 19b 
using the exoplanetar y atmospheric modeling and retrieval 
meth od developed in IMadhusudhan fc Seagerl l|2009l . l20ll 
l201lf ). We model a plane-parallel atmosphere of WASP- 19b 
observed in thermal emission at secondary eclipse. The day- 
side spectrum of the planet is generated using line-by-line 
radiative transfer, with constraints of hydrostatic equilib- 
rium and global energy balance, and includes the dominant 
sources of infrared opacity expected in gaseous atmospheres 
at high temperature. Our sources of opacity include molec- 
ular absorption due to H 2 Q, CO, CH4, CQ 2 , NH 3 , TiO, and 
VP (iFreedman. Marlev fc Loddersj 120081 ; iRothman et all 
120051 ; iKarkoschka fc Tomaskd boidl and H 2 -H 2 co l lision- 



induced abso rption (|Borvsow. Jorgensen fc Zheng! 1 19971 ; 
lBorvsow|[2003 ). The concentrations of the species and the 
pressure-temperature ( P—T) profile constitute the fre e pa- 
rameters in the model (IMadhusudhan fc Se agcr 2009). We 



explore the parameter space of th e model using a Markov- 
chain Mon te Carlo scheme (see IMadhusudhan fc Seagerl 
l20ld . l201ll ). and constrain regions of parameter space con- 
sistent with the measured planet-to-star flux density ra- 
tios at different levels of fit. Our goal is to constrain the 
existence of a possible temperature inversion, the dayside- 
to-nightside redistribution efficiency, the concentrations of 
the different molecular sp ecies, and the C/O ratio (e.g. 
IMadhusudhan et al.l l2011ah in the dayside atmosphere of 
WASP-19b. In what follows, we discuss model solutions that 
explain the data within the 1-a observational uncertainties, 
as shown in Fig. [4] 



4.2 Temperature inversion 

The data indicate the lack of a strong temperature inver- 
sion in the dayside atmospheres of WASP-19b. The ob- 
servations and two models fitting the data are shown in 
Fig. [4] The lack of a temperature inversion in WASP-19b 
is evident from the data even without detailed modeling. 
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Firstly, the brightness temperature in the 4.5 um chan- 
nel, Tb,4.5, is lower than that in the 3.6 um channel, 7b, 3. 6 
(Table [3} . The presence of a temperature inversion is of- 
ten indicated by a Tb,4.s value considerably higher than 
the 21,3.6 value, due to strong CO emission and some H2O 
emiss i on in the 4.5 urn channel dBurrows. Budai fc Hubenvl 
120081 ; iFortnev et at] |200S| ; iMadhusudhan fc Seagerl l20ld ) 
Secondly, the brightness temperatures at 1.6 and 2.09 urn 
are larger than those in all four IRAC channels, which are at 
longer wavelengths. Since the bands at 1.6 and 2.09 um are 
windows in molecular opacity, they are expected to probe 
temperatures in deeper layers of the atmosphere compared 
to any of the IRAC channels. Therefore, the high temper- 
atures in the 1.6 and 2.09 urn bands compared to all the 
IRAC channels imply temperature decreasing outwards in 
the atmosphere, and hence the absence of a temperature 
inversion or, at most, the presence of one too weak to be 
detectable with the available data. Two model P—T pro- 
files without temperature inversions and the corresponding 
model spectra are shown in Fig. [4] All the IRAC data can 
be explained by molecular absorption in the atmosphere due 
to the temperature decreasing outwards. 

The lack of a temperature inversion in WASP-19b of- 
fers a new constraint on existing classification schemes of 
irradiated giant exoplanets. WASP-19b falls in the category 
of highly irradiated hot Jupiters which have been predicted 
to host temperature in versions due to TiO and VO, assum- 
ing solar abundances (jFortnev et al.l [2008^ - the so called 
'TiO/VO hypothesis'. However, our finding of a lack of a 
strong temperature inversion in WASP-19b implies that ei- 
ther TiO and VO are depleted or an entirely different pro- 
cess is at play. If the composition is oxygen-rich, the lack of 
a temperature inversion in WASP-19b can be explained if 
TiO and VO are depleted in the upper atmosphere due to 
gravitational set tling, which can be significant if th e vertical 
mixing is weak ijSpiegel. Silverio fc Burrowsll2009l ). On the 
other hand, if the composi tion is carbon-rich, TiO and VO 
would be naturally scarce jMadhusudhan et aljfeoilri ). 

iKnutson. Howard fc Isaacson! l|2010l ) posit that the high 
UV flux impinging on those planets orbiting chromospher- 
ically active stars could destroy the high-opacity, high- 
altitude compounds that would otherwise lead to temper- 
ature inversions. With log i? HK = -4.50 ± 0.03, WASP- 
19 has a similar activity level to that of the handful of 
other stars around which hot Jupiters without tempera- 
ture inversions (logi?HK = —4.9 to —4.5) are known to 
orbit. Those planets thought to have inversions orbit qui- 
eter stars, with log -Rhk = —5.4 to —4.9. Thus, our find- 
ing that WASP- 19b has no inversion support s the hy- 
pothesis of IKnutson. Howard fc Isaacsonl (|2010T ) and use- 
fully adds to the handful of such systems known. Spitzer 
routinely measures the thermal emission of planets at 
3.6 and 4.5 um. For planets with temperature inversions, 
CO and water switch from absorption to emission, re- 
sulting in a higher flux in the 4.5 urn band, in which 
these molecules have features. IKnutson. Howard fc Isaacsonl 
(2010) proposed a model-independent, empirical metric for 
classifying hot Jupiters, which we denote with This is 
the gradient of the measurements at 3.6 and 4.5 urn, i.e. 
(AF4.5 — AF4.5)/0.9um, minus the gradient of the black 
body that is the best-fit to the two measurements. A positive 
(Rvalue would suggest an inverted atmosphere and a strongly 



negative C~ va hre would indica t e a no n-inverted atmosphere; 
IKnutson. Howard fc Isaacsonl (|2010T ) suggest a delineation 
around £ = —0.05 per cent um -1 . For WASP-19b we calcu- 
lated £ = —0.031 ±0.043 per cent um -1 . This value is at the 
border between inversion and no-inversion and therefore is 
consistent with our finding that WASP- 19b does not have a 
strong inversion. 

We note that in their activity-in v ersion plot (their 
figure 5), IKnutson. Howard fc Isaacsonl (|2010h omit XO-3 
but include TrES-4, HAT-P-7 and WASP-18, all of which 
have T ff > 6200. XO-3b has a temperature inversion and 
XO-3 has a \ogR' HK index indicative of activity, which 
seems to contradict the activity- inversion hypothesis. The 
other three planets have inv e rsions and orbit quiet stars. 
IKnutson. Howard fc Isaacsonl (|2010P ) concluded that in fact 
XO-3 is likely to be chromospherically quiet, based on a vi- 
sual inspection of their spectrum and having noted that the 
logR' HK calibration is unreliable for stars with T e g > 6200. 
Perhaps then TrES-4, HAT-P-7 and WASP-18 are also sus- 
pect since they also have T e g > 6200. 



4.3 Atmospheric composition 

We find that the observations can be explained by models 
with oxygen-rich as well as carbon-rich compositions. The 
absorption in the near-IR (1.6 and 2.09 um) bands is min- 
imal due to the lack of major molecular features. The con- 
straints on the composition come primarily from the IRAC 
data, which together encompass features of CO, H2O, CH4, 
and CO2. The near-IR data, however, are critical to con- 
straining the temperature of the lower atmosphere and thus 
are key in anchoring the model-atmosphere spectra to the 
measured SED. Two models with different C/O ratios, C/O 
= 0.5 (oxygen-rich) and C/O = 1 (carbon-rich), a r e show n 
in Fig. |4] As demonstrated in IMadhusudhan et al.l (|2011bT ). 
CO is a dominant carbon-bearing molecule in both C-rich 
and O-rich regimes. Consequently, the 4.5 um absorption in 
both models in Fig.|4]is caused primarily by CO absorption; 
in the O-rich model CO2 contributes additional absorption 
in this channel. The absorption in the 3.6, 5.8, and 8.0 mi- 
cron IRAC channels in the O-rich model is caused primarily 
by H2O absorption, whereas absorption in the C-rich model 
is caused by a combination of H2O and CH4 absorption; 
H2O is depleted by a factor of 100 and CH4 is enhanced by 
a factor of 1000 with respect to the O-rich model, both o f 
which are chemically feasible (|Madhusudhan et alj|2011b| y 
The principle difficulty in differentiating between the two 
models with the current data are the large uncertainties in 
the 5.8 and 8.0 um IRAC data. For example, a high 5.8 um 
point would indicate low water absorption, and henc e high 
C/O, as demonstrated in Madhus udhan et al.l (|2011al ). New 
observations in the near infrared can differentiate between 
spec tra from the carbon-rich an d oxygen-rich compositions 
(e.g. IMadhusudhan et al1l2011aT ). The water abundance can 
be measured via transmission spectoscopy of the 1.4 um 
water band using the G141 grism of HST/WFC3; t hese ob- 
servat ions were recently peformed for WASP- 19b bv lDemingl 
(2009). We could measure, or at least place useful constraints 
on, the TiO abundance with ground-based occultation ob- 
servations in the z and J bands. 
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Figure 4. Spectral energy distribution of WASP-19b relative to 
that of its host star. The black dots are our fitted planet-to- 
star flux density ratios. The observation-band transmission curves 
are shown along the abscissa. The planet-to-star flux density ra- 
tios indicate the absence of a strong temperature inversion. Two 
model-atmosphere spectra are shown, each of which fits the data 
and lacks a temperature inversion. One model has a C/O ratio of 
0.5 and is shown as a green line. The other has a C/O ratio of 
1.0 is depicted by a red line. The corresponding band- integrated 
model fluxes are shown as dots with the same colour scheme. 
The models' pressure-temperature profiles down to 100 bar are 
shown in the inset. Assuming zero albedo, the maximum dayside- 
to-nightside redistribution efficiency is 20 per cent for the C-rich 
model and 34 per cent for the O-rich model. The three dotted lines 
show planetary black bodies with temperatures of 1800, 2250 and 
2900 K (chosen to give a sense of scale to the spectra) divided by 
a stellar spectrum. 



4.4 Orbital eccentricity 

For a circular orbit, mid-occultation occurs half an orbital 
period after mid-transit. We find the occultation to occur 
only 4^30 s later than this and constrain both ecosoj and 
esinw to a small region around zero. Hence, the orbit is 
very nearly circular, though the available data do permit a 
small, non-zero eccentricity providing that the major axis of 
the orbit is near-aligned with our line of sight, such that the 
occulation time is not affected (i.e. |u;| ~ 90; Figure [SJ. We 
place a 3-er upper limit on eccentricity of e < 0.027. 

Hot Jupiters are considered to have been moved inwards 
to close orbits through planet-planet scatterin g or by the 
Koza i mechanism and tidal circularisation (e.g.: Naoz et al. 
20 111 ; IWu fe Lithwickll201ll ; iBatygin. Morbidelli fe Tsiganis 



2011f ). However, for the very shortest-period systems, such 



as WASP-19b, it is unlikely that they could have been 
moved directly to their current orbit, since that would 
have required careful fine-tuning to avoid destruction by 
collision with the star. Thus most likely WASP- 19b was 
first moved to an orbit near ~ 2 Roche radii and has 
since spiralled inwards through tida l orbital decay (see 
iGuillochon. Ramirez-Ruiz fc Linll201ll and t he discussion of 
WASP-19b specifically in lHellier et al]|201ll) . 

From tidal theory, the circularisation of a hot Jupiter's 
orbit is thought to proceed much faster than the infall, 
and this is consistent with the the observation that hot 
Jupiters tend to be in circular orbits. Thus the suggestion 



<a 0.02 




180 



/degrees 



Figure 5. The MCMC posterior distributions of e and ui. The 
white, blue and black contours are, respectively, the 1-, 2- and 3-cr 
confidence limits. The shading of each bin is proportional to the 
logarithm of the number of MCMC steps within. The available 
data exclude large values of e for any orbital orientation, and only 
very small values of e are permitted unless \lu\ 90. 



that WASP- 19b has undergone significant tidal decay, from 
~2 Roche radii to the current 1.2 Roche radii, leads to the 
expectation that the current eccentricity will be essentially 
zero, in line with our results. 
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